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We use entropy-energy arguments to assess the limitations on the running time 
and on the system size, as measured in qubits, of noisy macroscopic circuit-based 
quantum computers. 



Quantum information is essentially a nonequilibrium theory in the sense 
that its basic primitives are decomposed into temporal sequences of local op- 
erations that are performed on physical systems driven far from equilibrium. 
The requirement of locality becomes all the more important when we consider 
quantum information processing in macroscopic systems (multiqubit quantum 
computers), because local operations will be used to model both the compu- 
tation proper and the errors affecting it. 

An important issue to address is the stability of a macroscopic quantum 
computer (MQC) in the presence of noise. There are two notions of stabil- 
ity for quantum-mechanical statesti namely global and local thermodynamic 
stability. For a microscopic system (i.e., one with finitely many degrees of 
freedom) these are equivalent and amount to the following. Let H be the sys- 
tem Hamiltonian, (3 the inverse temperature. Then a state (density operator) 
p is thermodynamically stable (both globally and locally) if it minimizes the 
free energy functional -F/3(p) = tr {pH) — {1/ p)S{p), where S{p) = — tr (pin p) 
is the von Neumann entropy. However, the two stability notions are rather 
different for systems with infinitely many degrees of freedom. A state p is 
globally thermodynamically stable (GTS) if it minimizes the specific free en- 
ergy (i.e., free energy "per particle"), and locally thermodynamically stable 
(LTS) if no local modification of it yields a state with lower specific free en- 
ergy. It can be shown that any GTS state is also LTS, but the converse is 
generally false. 

States that are LTS but not GTS are referred to as metastable statesB 
An apt example comes from laser physics. Imagine an active medium con- 
sisting of a large number of three- level atoms with the states |0), |1), |2), each 
atom initially in the ground state |0). The atoms are pumped to the level |2), 
and then decay nonradiatively to This results in a relatively long-lived 
metastable state of population inversion. As far as MQC's are concerned, we 
are interested mainly in their metastable states, and the goal of error cor- 
rection (and of judicious MQC design in general) is to preserve metastability 
for the duration of the computation. Specifically, we can isolate two separate 
aspects of (meta) stability for MQC's — temporal and spatial. The former 

1 



refers to the maximum number of operations that can be carried out on a 
noisy MQC using given energy resources, while the latter is related to the 
maximum number of qubits that can be processed reliably on a iioisy MQC. 
We will address these two issues using entropy-energy arguments,a a standard 
technique in statistical physics. Due to lack of space, our presentation will 
be rather sketchy; we will supply the details in a separate paper.El We note 
that even though we consider macroscopic systems, there is no need to pass to 
the thermodynamic limit because we are concerned with local thermodynamic 
stability. 

First of all, let us fix a model of an A^-qubit quantum computer. We 
may imagine its operation as a sequence of alternating computation and noise 
steps. For a circuit-based QC, a computation step is an application of a tensor 
product of one- and two-qubit universal quantum gates, while a noise step is 
described by a general trace-preserving completely positive linear map T on 
density matrices. We will henceforth refer to such maps as channels (note 
that any unitary transformation is an instance of a channel). Specifically, if 
H ~ is a single-qubit Hilbert space, then the state of the A^-qubit computer 
is a density operator p on Ti,®^ . The composition of a computation step and 
a noise step is a mapping p TUJpU*), where U is unitary. We adopt 
the model of local stochastic noisejj so that T can be written in the form 
(1 — e) id+eR'^^ , where id is the identity channel, i? is a channel on 2 x 2 
density matrices, and e is a small positive number that quantifies the noise 
strength. 

We make two simplifying assumptions. The first one is needed for the 
analysis of the temporal stability and concerns the noisy channel R. Namely, 
we will take R to be bistochastic, i.e., R{1) = I, and strictly contractive Ja 
which means that there exists a constant C € [0, 1) (called the contrac- 
tion rate) such that, for any two density operators p,a, \\R{p) — R{a)\\-^ < 
C\\p— it||j^, where \\X\\-^ = tv\X\ is the trace norm. For instance, the much 
studied depolarizing channel, R{p) = (1 — X)p + AI/2 for some A G [0, 1), 
fits this description. Our second assumption, necessary only for the analysis 
of the spatial stability, states that the number of gates applied during any 
computation step is bounded by a positive number K that depends only on 
the particular algorithm, but not on the number of qubits n. This mmns that 
our analysis will be inapplicable to highly parallelized computation^ 

We will relj^Jipon the following two theorems, the proofs of which are 
given elsewhere.Bu 

Theorem 1. Consider the channel Rn — R"^^ , where i? is a bistochastic 
strictly contractive channel on 2 x 2 density matrices with contraction rate 
C. Then, for any 2^ x 2^ density matrix p and any positive integer m, we 
have S[R'}}{p)] - S{p) > ||p-2-^l||J, where S{p) = -tr(plnp) is 

the von Neumann entropy and II X II 2 = yJtT:{X*X) is the Hilbert- Schmidt 
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norm of the operator X. 

Theorem 2. Let pi, i — be n mutually commuting den- 

sity operators. Suppose that there exists a constant k > such 
that, for any i, J^j^ti^^ (PjPi) — Let p ~ '^"^X^ILiP*- Then 

Sip) > n-i ELi Sip^) + Inn - 2^. 

Let us first consider temporal stability. Suppose that our computer op- 
erates on N qubits, and that we are given "energy resources" E (this pa- 
rameter could be determined, e.g., from_the so-called "control Hamiltonian" 
representation of quantum computationEl) . We are interested in the entropy 
increase after m steps of noisy computation. Since entropy can only be pro- 
duced by the channel T and not by the unitarily implemented gates, the 
entropy gain is given by AS'(p, rn) = S[T"^{p)] — S{p), where p is the ini- 
tial state of the computer. Now T™ = J2T=o (7)(1 " e)'""^e^'i?jv, where 
Rn = i?^^, so that by concavity of the von Neumann entropy we have 
A5(p,m) > ^7^0 (7)(1 - e)'"-^eJ" {s[R%ip)] - Sip)}. Using Theorem 1 we 
can write 

AS{p,m) > Cp,JV E P (1 - e)"-^"e^"(l - C'^), (1) 

where (^p^N € [0, (2^ — l)/2^+^] depends only on p and N. Carrying out the 
summation in (|l|) we obtain 

ASip,m) > Cp,N {1 - [1 - e(l - C^)]'"} = Cp^meil - C^) + o(e), (2) 

where o(e) stands for terms that go to zero faster than e as e — > and therefore 
can be neglected when e is sufficiently small. The corresponding free-energy 
shift is AFf^ip, m) = E — (l//3)AS'(p, m), where (3 is the inverse temperature. 
It is now easy to see from the estimate (^) that AFp (p, to) will be negative for 
all TO > C,p^N,Tf3E, where Cp,n,t is a constant that depends on p, N, and T. In 
other words, the state of the MQC will not be metastable for large enough to, 
or, to put it bluntly, the MQC will "crash" after 0(/3i?) operations. However, 
because the produced entropy may be at least partially drained off by error 
correction, this estimate merely tells us how often it should be carried out. 

We pass now to spatial stability. We note that in a circuit-based computer 
each computational step is an application of a bounded number of universal 
quantum gates, and the initial state of the computer may be taken separable 
without loss of generality. Therefore at every instant of the computation the 
set {1, . . . , N} can be partitioned into L disjoint clusters Ci, . . . , Cl, so that 
the state of the computer has^the form p ~ (S'^Li p''''\ where /d^'-* is the state 
of the qubits in the Z'^ cluster .□ Let us suppose for simplicity that the clusters 
all have the same size d, so that the computer operates on iV = Ld qubits. 
Let us partition the set {Ci, . . . , Cl} into n disjoint blocks 5i, . . . , Sn, each of 
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which contains k clusters; thus N ~ nkd. Now, for each i ~ 1, . . . , 71, define Tj 
to be the channel that acts as a depolarizing channel a 1-^ (1 — A)cr + AI/2'^ on 
the state of each of the clusters in Si and leaves the remaining clusters intact. 
Then pi = Ti{p) are n mutually commuting density operators for which we 
have the elementary estimate 

trY,^p,p^<{n-l){l-\+^}j . (3) 

Note that for fixed n, A, and d we can always choose k so large that the 
r.h.s. of (||) is bounded above by unity. Thus the conditions of Theorem 2 are 
satisfied with k, — 1, and the entropy gain due to the channel T = n^^ ^^=1 Ti 
is O(lnn). Because the energy shift due to T depends only on the number K 
of gates applied during each computational step, the free-energy shift will be 
negative for large enough n. We may thus conclude that the number of qubits 
that can be processed reliably (that is, kept in a metastable state) in a MQC 
is bounded from above as 0{e^^), where /3 is the inverse temperature, and E 
depends on K. 

To summarize, we have shown that there exist upper bounds on the run- 
ning time and on the size (in qubits) of circuit-based noisy MQC's. However, 
these constraints become significant only when PE is small (as in, e.g., en- 
semble quantum computation using NMR). Furthermore, our analysis applies 
only to circuit-based computers without parallelization; massively parallel 
non-circuit rondels, such as the "one-way quantum computer" of Briegel and 
Raussendorf,cl are likely to be intrinsically thermodynamically stable. 
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